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<N ■ ABSTRACT 

The double meson production in ultraperipheral heavy ions collisions is addressed, focusing on 
q\ ■ the particular case of pJ/^ from two-photon reactions. The cross section at photon level is 
obtained using distinct parameterizations for the gluon distribution on the light meson. The 
^ ■ resulting estimates for the nuclear case are presented and discussed. As a by product, we esti- 
mate the double p production cross section using the Pomeron-exchange factorization relations. 

m 
o 

q ■ Relativistic heavy-ion collisions are a potentially prolific source of 77 collisions at high 

energy colliders [1, 2]. The advantage of using heavy ions is that the cross sections varies as 
Q_i! Z 4 a 4 rather just as a 4 as in the e + e~ collisions. Moreover, the maximum 77 collision energy, 
W^, is 2'j/Ra, about 6 GeV at RHIC and 200 GeV at LHC, where Ra is the nuclear radius 
and 7 is the center-of-mass system (cms) Lorentz factor of each ion. In particular, the LHC 
will have a significant energy and luminosity reaching beyond LEP2, and could be a bridge 
to 77 collisions at a future e + e~ linear collider. For the large values of energy at LHC the 
hard QCD Pomeron is presumably the dominant mechanism of production. Consequently, the 
estimate of the cross sections should to consider the QCD dynamics effects. Recently, we have 
proposed in Refs. [3, 4] to investigate QCD Pomeron effects in photon-photon scattering at 
ultraperipheral heavy ion collisions. In particular, in Ref. [3] we have analyzed the diffractive 
double J/ 1 ® production in 77 collisions, with the photons coming from the Weizsacker - Williams 
spectrum of the nuclei. For that process our results have indicated that future experimental 
analyzes can be useful to discriminate the QCD dynamics at high energies. Here, we extend 
the analysis of double meson production for the pJ/^ case using a perturbative approach and 
provide reliable estimates for the cross sections concerning that reaction. Since the Pomeron- 
exchange factorization theorem [5, 6] allow us to obtain the pp cross section in terms of the 
pJ/^! and J/^J/ty cross sections, we use our results for these cross sections to estimate the 
double p production in ultraperipheral heavy ion collisions. Moreover, since the double meson 
production can also occur in photo-nucleus reactions if the multiple interactions are considered, 
we compare our predictions for the two-photon process with that presented in Ref. [7] for the 
photonuclear case. 
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Let us start considering the pJ/^ production process at the photon level, that is 77 — > 
pJ/^f, with almost real photons. In the calculations presented here one follows the pioneering 
work in Ref. [8]. The goal of the present work is also twofold: investigate a discrimination 
among models for the gluon distribution on the light meson in an enhanced nuclear cross section 
and providing reliable estimates for the ultraperipheral nuclear cross section. The differential 
cross section is estimated in a similar way as the elastic J/\I/ photoproduction off the proton 
[9]. Therefore, it reads as [8], 

(W^t) =C a em g p — — '— 5 [xG p {x,M J/iSI /A)] explB pJ/ ^t) , (1) 

(11 O CXem "1 j /v]> 

where C denotes factors of corrections discussed later on. The two-photon cms energy is denoted 
by W 77 , where x = Mj^/W^ and Mj/q, is the heavy meson mass. In the small-t approxima- 
tion, the slope is estimated to be B pJ /^ = 5.5 ± 1.0 GeV -2 (See further discussion and Ref. 
[8]). The light meson-photon coupling is denoted by g 2 p = 0.454 and the heavy meson decay 
width into a lepton pair is Y ee . 

The process above was proposed as a probe of the gluon distribution on the meson xG rho 
and, by consequence, a constrain for the photon structure. The enhancement in the sensitivity 
by taking the square of those distributions in the total cross section could discriminate them 
in measurements at the future photon colliders. In Fig. 1, we present the distributions for 
the parameterizations considered here as a function of two-photon cms energy. They are the 
following: the LO (dashed line) and NLO (solid line) GRS [10] parameterization; the SaSID 
parameterization [11] (long-dashed line) and the phenomenological Regge motivated ansatz 
(dot-dashed line) proposed at Ref. [8]. The latter is given by, 

xG rho (x, M^/4) = x G rho {x , Mj /9 /4) )^ , (2) 

with xq = 0.1 and the Pomeron intercept is considered as up = 1.25, in agreement with the 
effective power in the HERA data. In further analysis for the nuclear case we allow for a 
higher intercept, in order to simulate a BFKL-like behavior, motivated by the studies in the 
double J/* production [3]. The normalization, x G rho (x , Mj^/4), is given by the NLO-GRS 
parameterization. The differences among the parameterizations are sizeable as a consequence 
of distinct effective exponent A. The GRS LO one presents the steeper behavior, followed 
by SaSID (A = 0.3038). The Regge motivated and GRS NLO parameterizations are more 
close since the intercept for the Regge ansatz is similar to the effective power of GRS NLO 
(A ~ 0.228). These behaviors will produce distinct results for the total cross section at photon 
level and in the nuclear case, allowing discrimination in future colliders. 

Let us now calculate the total cross section at photon level considering the models for 
the gluon distribution on the light meson referred above. In order to do so, we address the 
correction factor C. It accounts for several improvements to the process as partons skewness 
(off-diagonal) [12], QCD NLO corrections to the 7J/ 1 ]/ impact factor [13] and real part of the 
amplitude. Following [8], one has C = C^g Cnlo Crp where the correction factors are taken as, 

Cos ^ 1-2, Cnlo ^ H 7T , (3) 
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r ~ i.u 7r2A x dlog[:rG(x,Q 2 )] ^ 
4 alog(l/x) 

In Fig. 2 we show our predictions for the total cross section <Jtot(ll — ► P J/^) as a function 
of two-photon cms energy, extrapolated up to 1 TeV. The behavior on energy of the total 
cross section is strongly dependent of the choice of the gluon content of the light meson. In 
particular, a strong enhancement of the cross section is predicted if we assume the LO GRS 
parameterization, Cnlo = 1 an d disregard the real part of amplitude. The other predictions 
are similar to the results obtained in [8]. It is worth mentioning that this process has been 
also calculated in Ref. [14], considering a two-Pomeron model supplemented by the dipole- 
dipole picture. There, the result is similar to the Regge ansatz and lower than the other 
parameterizations considered here. 

Recently, the FELIX Collaboration has proposed the construction of a full acceptance de- 
tector for the LHC, with a primary proposal providing comprehensive observation of a very 
broad range of strong-interaction processes. In particular, studies of two photon physics in 
AA collisions has been discussed in its proposal [15]. There, the differential cross-section 
g?ct(77 — > ViV 2 )/dt, with V* = p, J/^f was parameterized in the form: 



da /W\ CviV * ( W\ 

-(ry^VM)=An( w ) exp(tB VlV2 + Ua' VlV2 ln—J , (5) 



where A VlV2 = 1.1 nb/GeV~ 2 , B VlV2 = 2.5 GeV -2 , C Vl v 2 = 0.8 and a' Vl y 2 = for the pJ/^ 
case. In Fig. 3 we compare the FELIX prediction (long-dashed line) with our results obtained 
considering the Regge motivated and GRS(LO) gluon parameterizations as input. Moreover, 
since in this paper we are considering -B p j/* = 5.5 GeV~ 2 , we also present the FELIX prediction 
for the cross section if we assume this value for the slope. We see that there is a reasonable 
agreement between the distinct predictions for the cross section at the photon level. Our 
conclusion at this level is that the forthcoming photon colliders could experimentally to check 
our predictions, which will shed light on the vector meson and photon structure. 

One comment is in order here. Our value for the slope is corroborated by the recent results 
for p and J/\& production at HERA and the proof of the factorization hypothesis for the slope 
parameters [6], which predicts that 



B w ^ pp (W) x B rp ^j/^ P (W) 
B PP (W) 



where B 1P ^ PP , _B 77 _^j/^ p and B pp are the slope of the diffractive peak for the reactions 7p — > pp, 
7P — > J/^p and pp — > pp, respectively. Data on 7p — > Vp in the energy range 50 GeV < W < 
lOOGeV gives B 1P ^ PP 11 GeV~ 2 and B 1P ^j/^ p w 5 GeV -2 . From pp reactions one obtains 
B pp = 10 — 12 GeV -2 . Consequently, assuming the validity of the factorization we obtain that 
B^^pj/^ ~ 5.5 GeV -2 is predicted, which is in agreement with the ansatz proposed in [8]. 

We are now ready to compute the double meson production in ultraperipheral heavy ions 
collisions. For two-photon reactions, the cross section for the process AA — > AA p J/^> will be 
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written as, 



f dC 

vaa^aa p j/v{s) = J dr (r) (T 77 _ p J/* (s), (7) 

where t = s/ s, s = is the square of the cms energy of the two photons and s of the ion-ion 
system, dC ^ T ^ j s the two-photon luminosity and o^^p j/^(s) is the cross section of the 77 
interaction (For details related to the numerical expressions for the effective photon luminosity, 
see Ref. [3]). Our approach excludes possible final state interactions of the produced particles 
with the colliding particles, allowing reliable calculations of ultraperipheral heavy ion collisions. 

In Fig. 4 are presented the results for the distinct parameterizations discussed above, 
considering A = Pb. For sake of completeness, we also include the Regge ansatz with a higher 
Pomeron intercept, up = 1.4, in lines of the studies on double J/ty production [3]. A upper 
bound is given by the SaSID parameterization and the lower one by the Regge ansatz with the 
lower Pomeron intercept uojp = 1.25. At the LHC energies y^s = 5500 GeV, the cross section 
takes values between 680 nb and 2.6 mb, showing the process can be used to discriminate among 
parameterizations. It should be noted that the ultraperipheral cross section is dominated by 
not so high two-photons energies. Although the deviations among the models are quite sizeable 
at very high energies, as shown in Fig. 2, in the nuclear cross section such deviations are less 
sizeable. This is directly related to the effective two-photon luminosity, which peaks at smaller 
Wy-f. Therefore, correct estimates should include careful treatment of the low energy threshold 
effects. The nuclear cross section is enhanced in relation to the e + e~ case, which provides values 
of hundreds of pb's in contrast with units of mb in the peripheral nuclear case. 

In Fig. 5 the results considering the FELIX ansatz for the two-photon cross section and 
different values of the slope are presented. For comparison the predictions obtained with our 
approach and the GRS(LO) input are also shown. We see that the large difference between 
the GRS(LO) and FELIX results for low energies in the two-photon case (See Fig. 3) implies 
that the values predicted by FELIX collaboration are ever larger than obtained using our 
approach, independently of the value of slope used. As a cross check, we also present in Fig. 5 
the predictions for pj /\& production in ultraperipheral heavy ion collisions obtained using the 
FELIX ansatz and the kinematical cut proposed in Ref. [15] . We have that our results agree 
with those presented in that reference. 

As a by product we can use our results for double J/^f production obtained in Ref. [3] and 
that presented above for pJ/^ production to estimate the double p production in ultraperipheral 
heavy ion collisions, using the factorization theorem to relate the distinct cross sections. Given 
the assumption that single Pomeron exchange dominates, the following relation among the 
distinct differential cross sections is predicted 



^ (77 _ VlVl) = I|ZL^ 



1 2 

(8) 



The usual approximation in the small t region, da/dt = A exp (£?£), where A = (da/dt) t=0 
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gives 



ah, - op) - {BpJ ^ 2 x Mil - p^/^)] 2 , 9) 

a{11 PP) ~ B pp x B mm a( 77 - J/* J/*) " (9) 

The slope of the diffractive peak for the 77 — > J/tyJ/ty process was estimated in Ref. [3] as 
being Bj/^j/^ = I /ml, where m c is the charm quark mass. Using the factorization theorem 
for the slopes, similarly to made in Eq. (6), we estimate B pp ks 12.1 GeV -2 . In the Tables 
1 and 2 we present our predictions for the double p production in 77 and ultraperipheral AA 
collisions, considering different scenarios for the QCD dynamics. We can see that there is a 
large range of possible values for the cross section, which implies that future experimental data 
are essential to constraint the dynamics. It is important to salient that our predictions agree 
with those presented in Ref. [15]. 

In Ref. [7] was emphasized that due to the large cross sections for vector meson pro- 
duction in photon-Pomeron interactions, the probability of having multiple interactions in a 
single nucleus-nucleus encounter is non- negligible. Therefore, it is possible that in double pho- 
tonuclear interactions occur the production of two mesons, generating a similar signal of the 
two-photon processes discussed here. There, the authors have presented some estimates of the 
cross sections for productions of vector mesons pairs in lead beams at LHC neglecting possible 
correlations between the multiple interactions (See Table V of Ref. [7]). Consequently, it is 
important to compare the predictions for the production of meson pairs in two-photon and 
photon-Pomeron interactions. Our predictions indicate that the a {A A — > AApJ/^) < 10 pb 
and <j(AA — > AApp) < 1250 pb, independently of the scenario considered. These values are 
almost one order of magnitude below of the predictions presented in Ref. [7], which implies 
that the rates of photonuclear multiple interactions would be an important background for the 
two-photon reactions, also for double meson production. In principle, an analysis of the im- 
pact parameter dependence should allow to separate between the two classes of reactions, since 
two-photon interactions can occur at a significant distance from both nuclei, while a photonu- 
clear interaction must occur inside or very near a nucleus. We salient that the experimental 
separation between the two classes of processes is an important point, which deserves more 
studies. 

As a summary, in this paper we have estimated the pJ/^f and pp production in two-photon 
processes at ultraperipheral heavy ion collisions. We have verified that if these measurements 
are feasible, a discrimination among the current parameterizations for the gluonic content of 
the mesons is possible. As a consequence, the analysis allows to constraint the partonic content 
to the photon structure function. Using the Pomeron-exchange factorization relations, the rates 
for double meson production were calculated and contrasted with the current predictions in the 
literature. The values found are in agreement with previous estimates. Concerning possible 
background processes, it was verified that the multiple meson production in photonuclear reac- 
tions constitutes the main source. The experimental separation between such processes should 
be taken into account, in particular focusing on impact parameter dependence and in rapidity 
distribution. 
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SCENARIO 


c( / y / y — ► dp) (pb) 


BFKL(LO) + GRS(LO) 


5 x 10 3 


BFKL(LO) + GRS(NLO) 


800 


BFKL(LO) + Regge (u P = 1.25) 


50 


BFKL(MOD) + GRS(LO) 


5 x 10 5 


BFKL(MOD) + GRS(NLO) 


80 x 10 3 


BFKL(MOD) + Regge (oop = 1.25) 


5 x 10 3 



Table 1: The double p production cross sections in 77 processes (W — 100 GeV). The notation 
BFKL (LO/MOD) corresponds to the treatment for the double J/^f production in Ref. [3] and 
GRS (LO/NLO), Regge, correspond to the treatment for the p J/ty production presented here. 



SCENARIO 


a(AA -> AApp) (nb) 


BFKL(LO) + GRS(LO) 


25 x 10 3 


BFKL(LO) + GRS(NLO) 


4 x 10 3 


BFKL(LO) + Regge (up = 1.25) 


810 


BFKL(MOD) + GRS(LO) 


125 x 10 4 


BFKL(MOD) + GRS(NLO) 


20 x 10 4 


BFKL (MOD) + Regge (ujjp = 1.25) 


405 x 10 2 



Table 2: The double p production cross sections in ultraperipheral heavy ion collisions at LHC 
(y/s = 5500 GeV) for PbPb. Same notation as in the previous table. 
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Figure 1: The gluon distribution on the light meson for different parameterizations. The solid 
and dashed lines are NLO and LO GRS, respectively. The long dashed line denotes the SaSID 
parameterization and the dot-dashed one the Regge motivated ansatz. 
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Figure 2: The total cross section for the process 77 — > p J/ty as a function of two-photon the 
cms energy. The distinct curves correspond to the different parameterizations. 
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Figure 3: The same as Fig. 2 considering the FELIX anzatz for the cross section and different 
values of the slope. For comparison our predictions using the Regge motivated input is also 
presented. 
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Figure 4: The total cross section for p J/ty production in two-photon ultraperipheral heavy ion 
collisions as a function of the ion cms energy (A = 208). The distinct curves correspond to 
the different parameterizations. 
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Figure 5: The same as Fig. 4 considering the FELIX ansatz for the cross section with (solid 
curve) and without (dot-dashed curve) the kinematical cut. For comparison our prediction 
using GRS(LO) input is also presented. 
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